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Abstract 

Let Q he a domain in a smooth complete Finsler manifold, and let G be the 
largest open subset of Q such that for every a; in G there is a unique closest point 
from dQ to x (measured in the Finsler metric). We prove that the distance function 
from dn is in C,^'"(G U dQ), k > 2 and < a < 1, ii dil is in C^'". 

1 Introduction 



In [1] we studied the singular set of viscosity solutions of some Hamilton- Jacobi equations. 
This was reduced to the study of the singular set of the distance function to the boundary 
of a domain Q — for a Finsler metric. The singular set was defined as the complement of 
the following open set 

G := the largest open subset of Q such that for every x in G there is a , . 

unique closest point from d^l to x (measured in the Finsler metric). 
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In [1] we stated that if dQ is in C'^''", k >2 and < a < 1, then the distance function 
from the boundary belongs to C''^^'"{GUdQ). Recently Joel Spruck asked to sec the proof 
for a Riemannian metric and pointed out that the result would imply that the distance 
function would be in C^'°^(G U dO). In this paper we provide a proof of that in the Finsler 
case. This paper can be regarded as an addendum to [1]. 

We present two proofs of the C*^'" result. We use the notation, as in [1], 



for the length of a curve $,{t). </?(^;v) is homogeneous of degree one in v. For fixed ^, 
the level surface ^{$,]v) = 1 is smooth, closed, strictly convex, with positive principal 
curvatures. 

The first proof uses very little of [1] and is essentially self-contained. The second proof 
uses some structure from [1] , and may be of some interest to some readers. 

We actually prove a more general result here, involving conjugate points from the 
boundary. 

Definition. Conjugate Point. Consider a point y on dVL, and consider the geodesic 
^{y, s) from y going inside f2 "normal" to dVt (explained below) with s as arclength. The 
conjugate point to y is the first point x on the normal geodesic such that any point x" on 
the geodesic beyond x has, in any neighborhood of the geodesic, a shorter join from dVt to 
x" than our normal geodesic to x" . 

Normal. A geodesic V from a point y e dVt is "normal" to dVt if for x on F close to y, 
the geodesic is the shortest join from dQ. to x. 

To obtain the regularity in G we prove a slightly more general result which is local on 
dVL. Namely, suppose C is a neighborhood on dO, of a point y and that the normal geodesic 
r from y to a point X in 1] is the unique shortest join from C to X. If the conjugate point 
to y is beyond X then there is a neighborhood A oi X such that the distance from C to 
any point in A belongs to C*^'". See Theorem 1 below. 

We shall make use of special coordinates introduced in section 3 of [1] about a given 
normal geodesic F, going from a point y G dVt into Vt. In these coordinates y is the origin 
and the normal to dVt there and is the geodesic F. Furthermore, in these 

coordinates, tp has the following properties, see (4.1)-(4.6) in [1]. Here Greek letters a, (3 
range from 1 to n — 1, and Latin letters i, j range from 1 to n. 




(/?(te„; e, 
<p^j {ten; e, 



■n 



n 




(2) 
(3) 
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(py^{ten]en) = 0, ^^n(te„;e„) = 1, (4) 

f^j^k{ten;en) = 0, (5) 

(te„; e„) = 0, (6) 

(^^j^n(te„;e„) = 0. (7) 

In these coordinates for y e near the origin the geodesic ^{y, s) from y "normal" to 
dfl there satisfies 

^{y,0) = Viy) 

where V{y) is the unique vector-valued function on dfl satisfying (here vly) is the Euclidean 
interior unit normal to dfl at y) 

■ V{y) ■ v{y) > 

< ip{y;Viy)) = l (8) 
. V„(^(?/; V{y)) is parallel to u{y). 

Using these special coordinates, near the origin, dfl has the form 

y={x'J{x')), /(0')=0, V/(O0=O'. (9) 

We assume that / e C"='", A; > 2, < a < 1. 
The result we prove is 

Theorem 1 Assume that the conjugate point of the origin on the geodesic F — {ten} is 
beyond e„, and that there exists a neighborhood C ofO' on dfl such that {ten \ < t < 1} 
is the unique shortest geodesic from C to e„. Then there exist neighborhoods A of e.^ and 
A of 0' on dfl such that for any X in A there is a unique y E A and geodesic from y to 
X which is the shortest join from A to A. Furthermore, if d{X) is its length, then the 
Jacobian of the map X — > {d,y) is nonsingular at en, and d lies in Ci^^ in A, and y lies 
in Cj,;^''* in A. 



2 Second Variation 

Consider one parameter family of curves T{e,t) from A to t > 0, with r(0,t) = te„. 
We look at the second variation of its length /[r(e, •)]. For t small, it is clearly positive 
definite. The first i for which it fails to be strictly positive definite is the conjugate point. 
For a t = t + 5, S > 0, then the second variation of curves to t cannot be semipositive 
definite, and there would then be a shorter connection from A to ten near P. 
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The standard computation of second variation yields 
if 



= ^(r4=o) - /...,(0')Tf (0,0)rf(0,0). 

e=0 



Here J is the usual expression of the second variation if the bottom point were kept at the 
origin. Namely, 

Note that r" and do not occur in J. 

3 First proof of Theorem 1 
3.1 

Recalling (9) we shall denote the normal geodesic from y = {x' , f{x')) hy X = ^{x', s); this 
is a slight change of notation. The geodesic ^ and depend smoothly on s and their initial 
data, while the initial data depend C*^"^'" on x'. To prove the theorem, it suffices to show 
that the Jacobian of the mapping (x', s) to X at (0', 1) is nonsingular. It follows that d 
and y belong to C*'"^ ". Since Vxd — Xg, it follows that Vxd is in C''~^'°' and hence d is 
in C'^''" — as Spruck pointed out to us. 

We now prove the Jacobian is nonsingular. 

Write X = Since X^(0', 1) = (0', 1), the Jacobian of the mapping {x', s) to 

X at (0', 1) is simply 

Assume M is singular, without loss of generality we may suppose that 

x;(oM) = o'. (11) 

We construct a perturbation r(e, t) oiT — {ten \ <t <1} such that ({t) :— t^I^^q satisfies 

^[c] = /.„.,(oor(o)c^(o). (12) 

3.2 

Consider the geodesic ^{Sei, t) of length 1 starting at {Se^, /(5ei)), < 5 small and "nor- 
mal" to dO, there. Set 

at) = ^MSeut) . (13) 
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By (11), 

C(1) = 0. (14) 
We obtain an equation for C(t) by differentiating the geodesic equation 

with respect to 5, and setting S — 0. We find 

d 

V9^i^i(te„;e„)C^ = — ((/?^«^,(te„; e„)C^) . 
Here we have used property (5) of our special coordinates. By (7) and (6), 

(^ee^C^ = ^ (^(f^cyisC^) . (15) 

We have 

r(o) = (5?. (16) 

In addition, 

v.,m. (17) 



C(0) = ^e(5ei,0) 



6=0 00 



5=0 



By the last formula in (8) we have 

V.v? ((5ei, /(<5ei)); \/((5ei)) ■ (ei + /.,(5ei)e„) - 0. 

Differentiating in S and setting (5 = 0, we find, using properties of our special coordinates, 

^,1,, (0'; e„) (0') + (0') = 0. (18) 

Now we introduce the perturbation T(e, t) as follows 

r"(e,t) = eC"(i), 

T"(e,t) = te„ + (l-t)/(eei). 

The definition of is just to ensure that r(e, 0) lies on dfl. 
According to (10), 

Integrating the last expression by parts we find, using (15) , (16), (17) and (18), 



□ 



= -^„i,.(0';e„)yi^(0') = /.,.,(0') = /.„.,(0')r,"(0,0)rf(0,0). 
It follows from Section 2 that the second variation is zero. 
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4 Second proof of Theorem 1 

Second proof of Theorem 1. For X near e„ and for small a' = (cti, • • • , cr„_i) e R^~^, 
let r = T{a',X) be defined by (p{X; {a',T)) = 1 and r(0',e„) = 1. Since (/9^,n(e„; e„) = 1, 
by the implicit function theorem, r exists as a smooth function of {a', X) near (0', e„). 
Let, as on page 111 of [1], = 'r]{a', X, t) be the unique smooth function of, with ip = ip"^, 

^eiv;v) ^ -^Aiiv^v), t<i, 

satisfying 

via', X, 1) = X, f,{a', X, 1) = {a', r{a', X)). 

As explained in the last two lines of page 108 in [1], r]{a',X,t) is a geodesic with t the 
arclength. 

In the special coordinates described in Section 1, dQ has the form (9) near the origin 
with / G C'^'", k > 2, < a < 1. Since {te„ | < t < 1} is the unique shortest geodesic 
from C to e„, we know that for X close to e„, there exists x' close to 0' such that the 
"normal geodesic" starting from {x' , f{x')) will reach X as a shortest join from C to X. It 
follows that for some a' close to 0' and t close to 0, we have 

/ r^(a',X,t)-(x',/(x')) = 0, 

\ f)^'{a',X,t) -V^'ix') = 0, ^ ^ 

where V{x') := V{x', f{x')) satisfies (8). Note that 1 — t is the distance from C to X. 
To prove Theorem 1, we only need to show that the left hand side of (19), denoted as 

LHS, has nonsingular Jacobian — — — at {t,a',x',X) — (0, 0', 0', e„). Indeed, this 

would allow the use of the implicit function theorem to show that for X close to e„ 
and in a neighborhood of (0,0', 0'), there exists a unique C*^-!'" solution (t, a' , x') = 
(t{X) , a' {X) , x' {X)) of (19). Thus, Theorem 1 follows as explained at the beginning of 
Section 3. 
Clearly, 

^<^*^)(o,o',o',.„) = ^ -iC''-") ^ 



a (2n — 1) X 1 column vector, 

9iLHS) (r7.„(0',e.,0))\_ 



1 

\0' 



( (C(0',e„,0)) 

(r);'JO',e„,0)) |, (20) 
^ (C(0',e„,0)) 
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a {2n — 1) X (n — 1) matrix, 

-/ 



-vy'(o') 



^^^(0,0',0',e„)= I |, (21) 



a {2n — l) X (n — 1) matrix, where / is the (n — 1) x (n — 1) identity matrix and W :— (Vj^^). 
Thus 

..<«,o,o^o^.,).Mr-v.(|o;;:-)_-^,,), 

where := «) and ^ := (ry^J. 
By the last hne in (8), 

i.e. 

^A{x', fix')); V{x')) + vAi^', fi^')); V{x'))f,^{x') = 0. 

Differentiating in Xa and setting x' — 0' we find, using properties of our special coordinates 
(4), (5) and (6), 

Dl,ip{0'; e„) • VV'{0') + D^f(0') = 0, (23) 

where D^,</? := ((/^^/j^m). 

We now evaluate |^(0',e„,0) and |^(0',e„,0). It is proved in section 4.4 of [1] that 
there exists a C^'^ function / near 0' satisfying 

m = 0, v/(o') - 0', 

(D'fiO') - D'fiO')) > 0, (24) 

r]{a',en,0)^{y',f{y')), (25) 

ria',en,0) = V''iy'), (26) 

where V{y') := V {{y' , f {y'))) is determined by (8) with / replaced by /, and y' = y'{(j') 
satisfies 

det (l^(O')) ^ 0. (27) 

Note that (27) is given by (4.9) in [1], while (24) follows from corollary 4.15 in [1] together 
with the fact that e„ is not a conjugate point. 



Differentiating (26) in (Tq and setting a' = 0' we find 

<(O',e„,O) = l>(O0|^(O'), 

i.e. 

g;(o',e„,o)=vnoo|;(oo. 

Differentiating (25) in ctq and setting a' — 0' we find 

^(0'.e„.0)^|^(0'). 

Since 

V wv'{o')^{o') -vr(o') ) \ VV'{0') -VV'{0') J \ I J 

we liave, by putting (28) and (29) into (22), 

Tfie proof of (23), applied to / instead of /, yields 

Dl,^{0'- en)VV'{0') + D'f{0') = 0. 

Thus, by (23) and (31), 

{-D'fm D^m) = ( (^(0'; e„) Vr (0') -^^,^(0'; e„)Vr(0') 

( Dl,ip{0';en) ) ( Vy'(O') -VV'{0') ) 

and therefore 

det (l)V(O') - /^V~(0')) = det Z),V(0'; e.) det ^^(^^^^ _vr (0') ) ' 
Since D^,<^(0';e„) is positive definite, we deduce from (30), (27) and (32) that 
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